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Abstract. We review different aspects of subaqueous sediment transport. We discuss the static threshold 
and its dependency with longitudinal and transverse slopes, as well as cohesion. We describe the different 
regimes of transport: erosion and momentum limited bed load, and suspended load. In all these cases, 
we derive the expressions of the saturation flux q sat and the saturation length L sa t and discuss their 
dependencies. 
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1 Introduction: the saturation length 
paradigm 

Whenever a sedimentary bed is sheared by a water flow 
of sufficient strength, the particles at the surface are en- 
trained into motion. The loss of static equilibrium of a 
grain at the surface of the bed occurs for a value of the 
driving force F proportional to the submerged weight of 
the grain: F oc (p s — Pf)gd 3 , where d is the grain diameter, 
p a is the mass density of the sand grain, and pf is that 
of the fluid. As F is proportional to the shear force rd 2 
exerted by the fluid (r = PfU 2 is the basal shear stress), 
the non-dimensional parameter controlling the onset of 
motion is the Shields number (9, which characterises the 
ratio of the driving shear stress to the normal stress: 



When the bed is composed by non-cohesive large grains, 
phcnomcnological transport laws relate the sediment flux 
q - i.e. the volume of grains crossing a unit transverse 
length per unit time - to the basal shear stress r. When 
the bed is composed of small cohesive or cemented parti- 
cles, one rather uses erosion laws, which relate the erosion 
rate ip — i.e. the volume of grains leaving a unit bed sur- 
face per unit time - to r. These two different cases can 
be described into the same general framework, introduc- 
ing the concept of saturation length. A relation between 
the flux and the shear stress is possible if an equilibrium 
state is reached, e.g. in a steady and homogeneous situa- 
tion. Hereafter we call g sat the equilibrium or 'saturated' 
sediment flux. g sat is then a function of r. Any change of 
the basal shear stress, for instance due to a non flat bed 
profile, leads to an actual flux which differs from g sat , and 
which generically relaxes towards equilibrium over some 
distance that we call L sa t. Then, a situation for which the 



actual flux is given by a transport law corresponds to the 
limit L sa t — > 0, whereas a saturation length larger than 
the system size corresponds to a purely erosional process. 
For many pratical situations - e.g. formation of ripples 
and dunes - a first order relaxation of the type 

L sat d x q = <7sat - q (2) 

is a good approximation of the full dynamics. 

In this paper, we review different aspects of subaque- 
ous sediment transport: (i) the static threshold and (ii) the 
different regimes of transport. Our goal is to make clear 
the mechanisms responsible for these transport properties 
in order to gain understanding beyond the phenomenol- 
ogy. More precisely, the paper is constructed as follows. In 
the next section we compute the threshold Shields number 
(9th from a miscroscopic and a continuum point of view. 
We discuss its dependency with longitudinal and trans- 
verse slopes, as well as cohesion. In section [3] we describe 
the regimes in which transport is limited by the erosion 
process of the bed or momentum extraction by the grains 
from the fluid. In section [4] we discuss the case of sus- 
pended transport. In all these cases, we derive the expres- 
sions of the saturation flux g sat and the saturation length 
i sa t and discuss their dependencies. 

2 Static threshold 
2.1 Discrete description 

The bed load threshold is directly related to the fact that 
the grains are trapped in the potential wells created be- 
tween neighbouring grains at the sand bed surface. To 
obtain the scaling law of the threshold shear stress, the 
simplest geometry to consider is a single spherical grain 
jammed between the two fixed grains below it. We first 
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consider the situation in which the cohesive forces between 
the grains are negligible and the friction at the contacts 
is sufficient to prevent sliding. 

The generic situation is a uniform shear stress PfU 2 , for 
which the velocity vertical profile is linear in the viscous 
regime: 



9 

uiz 



(3) 



where v is the kinematic viscosity of the fluid, and loga- 
rithmic in the turbulent regime: 



In 1 



z 
rd 
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where k is the von Karman constant and r the ratio of the 
hydrodynamic roughness to the grain diameter (zq — rd). 
The rescaled bed roughness r is on the order of r = 1/10. 
We introduce the effective velocity U of the flow around 
the sand grains and define for convenience its normalised 
counterpart: 

U = = = • (5) 

V(Ps/Pf - l )9 d 

For the sake of simplicity we will take for U the fluid ve- 
locity at the altitude z = d/2. Summing up viscous and 
turbulent contributions, a good approximation of the rela- 
tion between the shear velocity « 4 and the typical velocity 
around the grains U is given by 



In (1 + l/2r) 



-TP 



(6) 



At the threshold, the horizontal force balance on a grain 
of the bed reads 



^P(Ps 



2 j2 



Ps)9d i = -C dPf U^d 



(7) 



where p is the avalanche slope for sand grains. Cd is the 
drag coefficient , which is a function of the grain Reynolds 
number TZUd/v. In the viscous regime, the grain drag co- 
efficient Cd is inversely proportional to the grain Reynolds 
number (Cd — s 2 7?. _1 ); it is constant (C<j = C^) in the 
turbulent regime. Note that the turbulent drag is not only 
due to the fluctuations induced by the grain itself, as in the 
case of a grain falling in a fluid at rest, but also to those 
present in the turbulent flow. This situation has never 
been studied properly so far. Note also that the drag force 
is a priori different in the vicinity of the ground and that 
the lift force could be non-negligible when the grain is 
trapped at the surface of the sand bed. In order to avoid 
introducing two many parameters, we will not consider 
these effects here. With a good accuracy, the drag law for 
natural grains can be put under the form: 

C d = (CU 2 + sTl- 1 ' 2 ) 2 with K=™. (8) 

Typical values found in sedimentation experiments per- 
formed with natural sand grains are Coo — 1 and s ~ 5. 



At this stage, we introduce the viscous size d v , defined 
as the diameter of grains whose free fall Reynolds number 
UiaX\d/v is unity: 



d v = ( Ps / Pf - fl -VS. 



(9) 



It corresponds to a grain size at which viscous and gravity 
effects are of the same order of magnitude. From the three 
previous relations, we get the equation for U^h'- 



' j \ 3/4 , . v 1/2 



which solves into 
Z/ th = ' 
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(11) 



The corresponding expression of the static threshold Shields 
number is finally: 



9 th = 2 
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In the viscous regime, the above relations simplify into: 



In the turbulent regime, we get: 
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(14) 



As can be seen in figure [TJ the most interesting range 
of grain diameters, between 100 pvo. and 1 mm, is pre- 
cisely the zone of transition between the two regimes. The 
crossover grain diameter d co scales as: 



,4/3 



(mOo) 1 / 3 



(15) 



which is around 200 pm. 

The laminar-turbulent transition, indicated by the grey 
zone in figure [TJd, is determined by the Reynolds number 
and thus depends on the grain diameter d. In most of the 
situations relevant in geophysics, d lies in the cross-over 
between these two regimes. Figure [l] shows the compar- 
ison between experimental data for natural sand grains 
and the full model. There is no adjustable parameter in 
the model as one can independently measure the avalanche 
slope p — tan 32°, the drag coefficients and the sand bed 
roughness r = 0.1. 
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Fig. 1. (a) Threshold Shields number & as a function of the 
grain size d, for natural sand grains in water. Symbols: mea- 
surement by [21] (o) and by 12 (■). Solid line: Comparison 



with the model (equations 11 and 12 1, (b) Threshold shear 



velocity u t h as a function of the grain size d, for natural sand 
grains. Symbols: measurement by [21] (o) and by [12] (■). Solid 
line: Comparison with the model. Dotted line: asymptotic be- 
haviour in the turbulent regime. Dashed line: asymptotic be- 
haviour in the viscous regime. Grey zone: zone in which the 
viscous sub-layer is larger than the bed roughness: for grains 
smaller than d = 100 fim, the bed is hydraulically smooth for 
any shear stress u* above the threshold Uth- 
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Fig. 2. Dependence of the threshold Shields number (9" h on the 
bed inclination a. Symbols: measurements by [TJ] for natural 
sand grains (■) and by [14] using glass beads (o) in a narrow 
channel. Solid line: expectation of the model for natural sand 
grains. Dotted line: approximation by cos a + sma/p, with 
H — tan 32°. 



2.2 Continuum description 

The transport threshold can alternatively be described us- 
ing a continuum description of two-phase flows. We con- 
sider here the viscous limit. The turbulent case can be 
obtained in a straightforward way. In a steady state, the 
equilibrium equation reads: 

p p <Pg + \7 -cr s + <PS7 -er f + -s 2 <P^U = (16) 

4 cr 



p f {l - $)g + (1 - <P)V • a f - - s 2 <P^U = (17) 



where <r s and cr? are respectively the sand borne and fluid 
borne stresses. Along the vertical direction, this equation 
reduces to the standard hydrostatic condition: 

P f = -Pf9 z et P s = -a s zz = -(/0p - Pf)$gz 
Along the horizontal direction, the force balance reads: 



dai 



(1-*) 



dz 
dot, 
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dz 4 d 2 



(18) 
(19) 



with a s xz = /icrf 2 and a{. z = r]U'(z) . Integrating the sec- 
ond equation, one gets the fluid velocity field: 



V 



rd /4(l-<2>) 
3s 2 <Z> 



exp 



/ 3s 2 z 
4(1 ~<P) d 



The fluid velocity relaxes exponentially to inside the 
sand bed, described as a porous medium. The first equa- 
tion gives the sand borne shear stress: 



1 — exp 



I 3s 2 $ z 
4(1 -<£) d 



As expected, a s xz tends in depth to the total shear stress 
t. The Coulomb criterion <J s xz j a s zz = \i is reached for the 
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first time in z = and gives a threshold Shields number 
equal to: 

/4#(1 - ¥) 
3 s 2 

For natural grains (s ~ 5), the result is quantitatively very 
close to that given by the discrete approach (0 t h — 0.07), 
and to the measurements in the viscous regime. 

2.3 Effect of longitudinal slope 

Due to gravity, it is much easier to transport sediments 
along the lee side of a sand bump than on the stoss slope. 
This effect is directly encoded into the dependence of the 
saturated flux on the transport threshold, which depends 
on the slope. Let us note a the bed angle at the scale of 
the saturation length. The force balance at threshold now 
reads: 

^(Ps ~ p f )gd 3 (sin a + p, cos a) = ^C d p f (U? h d) 2 . (20) 

We keep the notations u t h and (9th for the threshold shear 
velocity and the threshold Shields number for an horizon- 
tal bed (a = 0) and denote by u" h and 0" h the thresholds 
in the case of an inclined bed. 

The dependence of the threshold on the slope has been 
measured experimentally by |12j for sand grains, in the 
turbulent regime. The model matches quantitatively the 
results, without any adjustable parameter (solid line on 
figure [2]d) . Neglecting the fact that Cd is a function of the 
grain Reynolds number, the threshold shear stress can be 
written as: pj (u" h ) 2 = p/u 2 h (cos a + s'ma/p). This ap- 
proximation still gives a reasonable fit of the data (dotted 
line in figure |2|d). 

The slope dependence of the threshold shear stress 
has also been measured in the case of aeolian transport, 
for rough sand grains ([18]) and is again perfectly fitted 
by the model. However, there exists an unexplained dis- 
crepancy in the case of spherical glass beads in a viscous 
liquid (Fig. [T]d), as one would need a very low effective 
friction p ~ tan 70° much lower than the avalanche slope 
(p ~ tan 24°) to fit the data ([II]). With such a value, the 
gravity effect becomes completely negligible. 



At the threshold, the velocity V vanishes so that: 

g _ 4k 2 (pt + d x Z e x + d y Z e y ) 

~ the * " 3C d ln 2 (l + l/2r) y/1 + (d x Z) 2 + (d y Z) 2 ' 

(23) 

Defining the threshold Shields number for a flat bottom, 



= 4/iAt 2 

th 3C d ln 2 (l + l/2r) ' 



we obtain: 



yJl-(d y Z /pf + d x Z/p 
Ol ~ y/T+Jd x Z) 2 + (d y Z) 2 
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(24) 



(25) 



The threshold depends linearly on the longitudinal slope 
but quadratically on the transverse one. 



2.5 Effect of cohesion 

We now consider the case of cohesive granular media con- 
stituted of grains sufficiently small to lead to van der 
Waals interactions of the same order or larger than grav- 
ity. The adhesive force is related to the real surface of con- 
tact between grains. It depends on the loading history or, 
more precisely, on the maximum normal force exerted pre- 
viously. When the cohesive grains have gently sedimented 
at the surface, this force is simply the weight of a single 
grain. On the contrary, when the grains have stayed under 
a very high pressure or are surrounded by clay, the whole 
surface contributes to the adhesion force. Introducing the 
surface tension 7, the maximum adhesion force is of the 
order of tt 7 d, which leads to a modified threshold Shields 
number: 



£ ( 7 ) - 0th (0) 1 + 



67 



m(p« - Pf)gd 2 



(26) 



The grain size at which the cohesive effects arc of the same 
order as gravity is of the order of few millimeters in the 
case of capillary bridges and few tens of centimeters for 
clay. 



2.4 Effect of transverse slope 

We consider a solid slide (a grain) on a surface Z(x,y) 
inclined with a longitudinal slope d x Z and a transverse 
slope d y Z. In the turbulent regime, for which the grain 
drag coefficient is roughly constant (Cd — 1), the drag on 
a grain at the velocity Vt reads: 

^C dP fd 2 \Ue x -Vt\(Ue x -Vt), (21) 

and should balance the gravity and the solid friction: 

n (p s - pj)gd 3 



6 + (d x z) 2 + (d y zy 



(pt + d x Ze x + d y Ze y ) 



(22) 



3 Bed load, from erosion to momentum 
limited regime 

When a sand bed is submitted to a flow, only a small 
fraction of the grains at the surface are entrained and this 
erosion process takes some time to occur. Close to the 
threshold, there is a regime in which the concentration of 
mobile grains is not sufficient to induce a significant reduc- 
tion of the flow strength in the transport layer. The grains 
transported in this dilute situation are isolated from each 
other. The saturation of the flux is then controlled by the 
erosion rate and the disorder of the sand bed. This idea is 
due to [9], who proposed a semi-phenomenological model 
that uses extensively experimental results. As the concen- 
tration of transported grains becomes larger, the transport 
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becomes limited by the available momentum: each time 
the flow entrains a grain from the bed and accelerates it, 
the grain exerts in return a stress on the fluid (0). The 
fluid in the transport layer is thus in equilibrium between 
the driving shear stress PfU 2 , the fluid-borne basal shear 
stress pfUf and this sand-borne shear stress. We propose 
here a model that includes these two regimes. 



3.1 Trajectory of a single grain 

As for the computation of the threshold in the previous 
section, we shall start with the discrete description of the 
trajectory of a single grain over a reg ular bed (Figj3f. To 
describe the elementary jump by one grain diameter, we 
write the equation of motion of one grain dragged by a 
flow of effective velocity U and submitted to gravity. We 
assume that it looses completely its energy during colli- 
sions, due to the thin layer of fluid between the grains. 
Contrarily to [7] and [5], we do not introduce any 'effec- 
tive friction': instead, we take into account the geometry 
of the trapping grains (Fig. [3]), which is the physical effect 
responsible for this friction. 

Starting from the position of the grain dug, we get 
its velocity V = d9u g+7r /2 an d its tangential acceleration 



d9u 



6+TT/2 



We assume that the reaction of the substrate 



is normal to it. Using again das a unit size and [(p s /pf — 
Vjgd] 1 / 2 as a unit velocity, the dynamical equation reads: 



CU 2 [0 2 - 2 cos 6U 9 + U' 



1/4 



d„ 
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3/4 
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(27) 
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9 2 - 2cos9Ud + U' 



1/4 



-U, 
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where p, has been taken equal to tan(7r/6) — 1/V3. This 
equation can be integrated numerically between 9 = —tt/6 
and 9 = tt/6, starting from 9 = 0. Note that there is no 
adjustable parameter in this description. 

In order to get the asymptotic behavior analytically, 
we perform an expansion with respect to 9 around —tt/6 
and with respect to hi around U^. We get a linear equation 
of the form: 



1 



2n(ht-U th ) 



3 3 / 2 aU, 



th 
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(28) 
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Fig. 3. Geometry of the trajectory calculation. 



The solution of this linear equation is: 

TT TT (Uth - U) 



1 



6 SaUth 
where r_ and r + are the solutions of the equation: 

„2 , ( 1 



(30) 



ih 



The time T needed for a grain to move by one diameter 
thus scales around the threshold as: 



T 



1 



In 



th 



U-U th 



p f d 



\j (Ps - Pf)g 
( 31 ) 

This erosion time diverges at the threshold as: 



Tex - In [0- 0th]- 



(32) 



On the other hand, T should resume to d/U far from 
the threshold. Integrating numerically equation (27 1, we 
have found convenient approximations that follow the cor- 
rect asymptotic behaviors just above the threshold and far 
from it. In the turbulent regime, it reads: 



T cx In 



Pfd 



(ps - Pf)g ' 



(33) 



where the pref actor is around 1.5. In the viscous regime, 
it reads 



T cx In 



6> + <9 



th 



e-e 



th 



d„ 



3/2 



p f d 



(p s - Pf)g' 



(34) 



with a prefactor around 17. 



3.2 Saturated flux 



We define Af(0) the fraction of grains at the surface sus- 
ceptible to be entrained at a Shields number 0. In the 
absence of flow, for = 0, N is null. It increases very 
quickly around the mean threshold Shields number (9th 
and reaches 1 at a second threshold 0m- N(0) reflects 
and encodes the distribution of potential wells at the sand 
bed surface. A simple choice is: 



N=0 if 0<0 mi 

0-0r, 



Af 
N 



0M - 

1 if > M 



if m < < 



M ■ 



(35) 
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Fig. 4. Left: mean grain velocity V = {{ps/pf — l)<?d) _1,/2 4- 
as a function of the rescaled Shields number 0/Oth m the 
turbulent regime (large d, solid line). Right: same, but in the 
viscous regime (small d, dotted line), with an further factor 

m s/2 • 



As suggested by the experiments of [9], a granular bed 
prepared by sedimentation is initially very disordered and 
consequently m is well above 6>th an d large. It then takes 
a very long time for the surface to re-arrange, leading to a 
drift of the distribution M toward larger and larger thresh- 
old shear stresses. This long transient is probably not rele- 
vant to geophysical situations in which the time and length 
scales of the systems are always sufficiently large to ensure 
that an equilibrium state for the geometrical arrangement 
of surface grains has been reached. 

When one grain is entrained, we hypothesise that there 
can be no further erosion in a surrounding area d 2 un- 
til the grain has moved by a distance comparable to its 
own diameter d. The erosion time T previously computed 
leads to a mean velocity that goes to at the thresh- 
old like 1/ In [0 — 6> t h] ■ This corresponds to a bifurcation 
much sharper than the standard saddle node bifurcation in 
\J0 — Oth- The experimental observation ([UlEllH]) that 
the grain velocity - i.e. the mean velocity conditioned by 
the fact that the grain moves significantly at the time 
resolution of the instrument - is non zero just above the 
threshold is thus consistent with this expression. Far from 
threshold, the grain velocity is just proportional to the 
flow velocity ([TUT^llS]). The length of a trajectory is in- 
dependent of the trap in which the grain was initially at 
rest and fluctuates around a well defined average value C. 
The mean area explored by the grain is Cd and contains, 
by definition of C, a mean number of potential wells suffi- 
ciently deep to trap the grain equal to 1. We then obtain 
the relation: C (1 -Af(0)) = d. 

The particles are entrained by the flow at a shear ve- 
locity Uf reduced by the presence of other particles. The 
relevant Shields number should thus be based on Uf 
rather than it* , to take the negative feedback of the trans- 
port on the flow into account. We consider the regime in 
which the grains are transported at the surface of the bed 
i.e. do not form a surface sheet flow (|16J). They leave 
the bed with a velocity and collide back the sand bed 
with a velocity . The sand-borne shear stress is propor- 
tional to the sand flux and to the difference (i^ — Uf). The 
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Fig. 5. (a-c) Saturated flux g sat rescaled by y — lj gd 3 as 

a function of the Shields number 0: raw data (□) gathered by 
[13] : same after local averaging (•); raw data (o) obtained by 
[12] , The solid line is the best fit by the model described here 
(equations |38| and |40[ | , which gives: (9 t h = 0.045 and (9m = 
5<9th- (a) Lin-lin plot; (b) Log-log plot; (c) Log-log plot as a 
function of — 0th- (d) Velocity Uf in the transport layer as 
a function of determined from the model (equations 38 and 
40}. 



equilibrium of the fluid in the transport layer leads to the 
relation valid in the saturated state (q = q sa t) but also 
during the transient of saturation 



Pf u * = Pf u f + Ps4> — - r — 1- 



(36) 



The grain trajectory is encoded into a single quantity, 
(vi — vi)/C, which is a function of the reduced shear ve- 
locity Uf and not of . 

The saturated flux q sat is simply the erosion rate times 
the hop length C = d/ (1 — Af(0)). It can be expressed as 
an integral over all possible local configurations: 



<Zsat 



nd 2 

60(1 —Af(0)) j,. 



N' (gth) 
T(G th ,9) 



d0. 



(37) 



where <fi is the bed volume fraction. A good approximation 
of this integral is given by: 



<7sat oc 



(j) (u 2 



M 



U 2 ) 



In 




where c m ~ 1.41 in the turbulent regime and c m ~ 0.605 
in the viscous regime. In particular one can recover from 
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this expression the expected asymptotic behaviours. By 
construction, it vanishes at the shear velocity u t h and di- 
verges at um ■ For the sake of simplicity, we limit the dis- 
cussion to the turbulent case. The transposition to the 
viscous case is straightforward. Assuming that the veloc- 
ity increment V[ — U| scales on the grain mean velocity, the 
ratio (vi —v^)/C is proportional to the hop time T, which 
is also the erosion time. Equation (36 1 can be solved to 
get the flux: 



<7sat oc 




.(39) 



Eliminating the flux between the equations ( 38 1 and ( 39 1 
one obtains a relation between Uf and . 



(ul - u 2 h ) cx (u 2 - u 2 th ) 
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4> (u M - u 2 ) 



In 
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111 



Uf +"th 
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Just above the threshold, the flow in the transport 
layer is undisturbed: Uf ~ it*. This corresponds to the 
erosion limited regime and the flux can be approximated 
by: 



Given these reservations, the fit by the model derived here 
provides a good description of existing data. In particu- 
lar, the asymptotic behaviors close to the threshold and 
for large shear velocities are well captured. 



3.3 Saturation length 

The sand flux is the product of the grain velocity by the 
concentration of mobilised grains. Two important mech- 
anisms can thus control the saturation length: the length 
needed by the grain to reach its asymptotic velocity - the 
so-called drag length - and the length needed for erosion 
to take place i.e. the trajectory length C. The modelling 
of the drag length is a difficult problem as the trajectory 
takes place in a turbulent flow whose fluctuations are not 
due to the motion of the grain itself. The problem is thus 
very different from that of a sphere moving in a fluid at 
rest, a problem for which the drag law is calibrated. To the 
best of our knowledge, the motion of a sphere whose diam- 
eter lies in the inertial range of the turbulent flow is still 
)an open problem. We are thus left with the standard drag 
force formula ^CdP/U 2 d 2 , with a drag coefficient of order 
one. Then, solving the equation of motion, one obtains a 
drag length Ldrag around 2(p s / Pf)d. 

In summary, we consider here that L sa ± can be either 
limited by erosion, in which case it is expected to scale as: 



1 



<7sat OC 



hi 




l-Af(0) 



(43) 



Far above the threshold, the velocity profile inside the 
transport layer becomes independent of . The shear ve- 
locity Uf tends to mm, value for which all the grains at the 
surface can be mobilised. The flux then scales as: 



At small 0, L sa t is simply one grain diameter d\ at 0m, 
J\f tends to 1 so that L sat would diverge. The divergence 
of the trajectory length has been directly evidenced ex- 
perimentally in the viscous case, by [8]. However, the sat- 
uration length becomes limited by the grain inertia below 
0m, in which case we get: 



<7sat oc 



p f d 



Pf 



Ps V (Ps - Pf)a 



(42) 



2^d 
Pf 



(44) 



Thus, in the momentum limited regime, q scales at large 
velocities as u 2 and not ul as usually obtained. The same 
result is valid for aeolian transport, for which the nega- 
tive feedback of particles on the flow has been directly 
evidenced experimentally. Although many authors have 
followed Bagnolds, the scaling law in u 2 is, in that case, 
the best description of existing data [I] . In figure [5] we 
have re-plotted the measurements of the saturated flux 
collected by [T3] for the subaqueous case. It exhibits a large 
dispersion due to several factors. First, data coming from 
systems with different grain size distributions have been 
plotted together without any distinction of symbols. Sec- 
ond, the reproducibility of experiments is made difficult by 
the problem of granular bed preparation: we emphasise 
again that sedimentation leads to an out of equilibrium 
situation that can last for days. Third, it is difficult to 
estimate the basal shear stress in flume experiments due 
to lateral boundaries. We have also plotted the measure- 
ments performed by |12j . which are much less scatterred. 



One can see that these two predictions are difficult to test 
and discriminate. The erosion length gently increases with 
the shear velocity while the drag length is independent of 
it; the drag length increases with the density ratio pal Pf 
but, experimentally, this parameter cannot be easily var- 
ied by a large factor. 



4 Suspended load 

We consider a flow of thickness H above an erodiblc bed. 
We introduce </>, the mass concentration of particles in 
suspension. We limit the discussion to dilute turbulent 
suspensions, in which there is no significant feedback of 
the particles on transport. Assuming that the suspended 
sediment transport is linearly driven by the concentration 
gradient through a particle turbulent diffusion coefficient 
-D, the equation governing the evolution of <p reads 

d t 4> + u x d x ct> = d x {Dd x <t>) + d z {Dd z <t> + <t>V fal i) , (45) 



8 



B. Andreotti, P. Claudin: On the different regimes of subaqueous transport 



where Vfaii is the settling velocity, assumed uniform and 
equal to that in the dilute limit. Note that we neglect here 
the vertical component of the advection. The sediment 
flux q introduced above is here defined as 



q = — / on r <li . 

Ps 



(46) 



For simplicity, we consider a flow with velocity U assumed 
uniform in the vertical z-direction. We assume the friction 
velocity it* is assumed to be proportional to the flow ve- 
locity, i.e. U = Au*. The diffusion flux is assumed to be 
proportional to the gradient of the particle concentration, 
with diffusivity D = )Cu*H, where K. is a constant. 

In the saturated state, steady and homogeneous in the 
longitudinal direction, the vertical diffusive flux of parti- 
cles is balanced by the settling flux. With the boundary 
condition of zero flux at the free surface, equation (451 



governing the evolution of the concentration then reduces 
to 



D 



dz 



■Vfwi = 0, 



(47) 



so that the equilibrium particle density, denoted <P, de- 
creases exponentially with height: 



<P = •Pq exp 



(-4) 



with 



Vfall 



(48) 



In a thin layer close to the bed, of thickness of a few 
particle diameters, the settling velocity vanishes and the 
particle concentration increases to that of the fixed bed, 
(j) M sa 0.6. Just above this layer, which may be consid- 
ered as the bedload layer, say at z = 0, continuity of the 
erosion flux gives the boundary condition 



dcf) 

ipM = —D— at z — 0, 
az 



(49) 



which determines the expression of the concentration <Pq: 



$n = 



f(j) 
Van ' 



(50) 



Finally, the saturated particle flux per unit width, nor- 
malized by the water flux UH, is given by 



<?sat 

UH 



1 

UH 



H 



<PUdz = 



1 - c" Q 



(51) 



with the assumption that particles have negligible inertia 
and follow the stream. For small a, this dimcnsionlcss flux 
tends to ^o, as expected. 

The saturation length is controlled in this regime by 
the length needed by a particle at a typical height ~ H 
above the bed to settle down: 
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